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Eo INDRERUCTION 


The ocean and atmosphere behave as a coupled system with anomalies 
in one medium inducing anomalies in the other. Bjerknes (1966, 1969, 
1972) and Namias (1959, 1969, 1972) have investigated this interaction by 
relating observed sea-surface temperature (SST) anomalies to anomalous 
behavior in the atmosphere. Other investigations have been conducted 
with numerical simulations to test the response of the atmosphere to a 
prescribed ocean (Rowntree, 1972), the response of the ocean tc a pre- 
Scribed atmosphere (Haney, 1974) or the simultaneous effects of a coupled 
system (Manabe, 1969; Bryan, 1969; Wetherald and Manabe, 1972). Through 
observations and numerical simulations a greater understanding of the 
air-sea system is sought with improved forecasting of the atmospheric and 
Seeanve Variability an ultimate goal. 

The numerical model used in this investigation is that of Haney (1974) 
with subsequent improvements involving seasonal atmospheric forcing, para- 
meterization of convective mixing and surface wind stirring (Haney anc 
Davies, 1976) and nonlinear eddy viscosity (Haney and Wright, 1975). The 
goal of this investigaticn is to further improve this model by introducing 
salinity as a new variable. The model, based on the hydrostatic and 
Boussinesg approximation, is a 10-level primitive equation model in a 
closed rectangular basin with uniform depth of 4 kilometers. The basin 
which extends from the equator to 65°N and spans 90° of longitude is an 
idealized representation of the North Pacific Ocean. The horizontal re- 
solution of approximately 2° of latitude and 3° of longitude is provided 


by 33 uniformly spaced grid points in the north-south and east-west 


IG, 





fieeectens. In the vertical, the ten levels are located at 10, 32.5, 
62.5, 102.5, 162.5, 462.5, 900, 1700 and 3100 meters below the surface. 
Wind and differential heating which drive the ocean circulation are de- 
rivea from prescribed but time varying climatological values of solar 
radiation, cloudiness, surface air temperature, surface specific humid- 
ity and winds. Until the present study, the density in the model was 
assumed to be a linear function of temperature only. 

The model has been developed and is being improved in order to study 
large-scale thermal anomalies observed in the Central North Pacific 
Ocean by the North Pacific Experiment (NORPAX). 

During a "spin-up" phase, the model is integrated for approximately 
240 years of simulated time to produce a "model climatology". The result- 
ing circulation and temperature pattems adequately represent the large- 
scale north-south temperature variations and the major ocean gyres; 
however, the Alaskan gyre is not depicted. Once the model climatology 
has been established, temperature anomalies are injected into the model 
which is then integrated forward in time with or without anomalous atmo- 
spheric forcing. The evolution and migration of these anomalies can he 
related to their horizontal or vertical extent and atmospheric conditions 
to determine the controlling processes. Investigations by Hunt (1975) 
and Shiver (1977) detail the use of this model in such studies of large- 
scale ocean variability and the major conclusions of these works are pub- 
lished in Haney, Shiver and Hunt (1978). 

The motivation for this investigation is to determine the effects 
of adding a prescribed salinity cross section to the model. The salinity 
cross section is steady state and independent of longitude. This new 


Salinity field modifies the density field through a newly introduced 
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equation of state. The resulting modifications to the density field 
affect the hydrostatic pressure, and hence the horizontal pressure gra- 
eeeent, a5 Well as the vertical stability. 

This introduction of a prescribed salinity cross section is expected 
to improve the model in several ways. First, as noted above, the present 
model climatology (without salinity) does not depict an Alaskan gyre. 
Part of this model's deficiency may be due to the somewhat idealized 
form of the wind forcing at high latitudes; however, temperature varia- 
euens are Small and salinity can play an important role in the density 
structure and hence, the geostropnic flow. Thus it is expected that 
Salinity may improve the model's circulation at high latitudes. The de- 
Sirability of adequately simulating the Alaskan gyre in this meteorologi- 
cally significant area is apparent. 

In his work using the model in which density depended on temperature 
Suey, aunt (1975) noted that cold SST anomalies could not be maintained 
in the winter at high latitudes but were immediately convectively mixed 
downward. Cold SST anomalies are a physically realizable condition dur- 
ticewtater in high latitudes. Surface cooling occurs but the increasin 
Salinity with depth provides stability and maintains the cold surface 
anomaly. However, free convection occurs and must be modeled whenever 
the surface layers become more dense than the underlying layers. The 
model previously determined stability solely from the temperature profile. 
If the temperature decreases with depth (positive profile), the profile 
is stable and convection does not occur. On the other hand, an unstable 
profile (produced for example by a cold surface anomaly) is instanta- 
neously adjusted to neutral thereby destroying any cold surface anomalies 


Wieeqmoreducea tie instability. In reality, density, and therefore 


We 











Stapility, 1s a complex function of temperature, salinity and depth and 

a negative temperature profile can be stabilized by a positive salinity 
profile. In this manner, cold SST anomalies can be formed and maintained 
and one of the goals of this study is to make the model more realistic 

in this regard. Similarly, but with less apparent significance, a warm 
anomaly can be maintained below the surface layers. The inclusion of 
Salinity in the model can play an important role in the study of large- 
scale thermal anomalies by improving the climatological circulation pat- 
tern and altering the vertical stability, and hence, the convective 


adjustment parameterization. 
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It. DATA COMPILATION AND HANDLING 


The salinity cross section used in the model is steady state and in- 
dependent of longitude. The steady state assumption was necessitated 
by the lack of data. A seasonal variation in surface salinity tied to 
the variation in evaporation is known to exist in several geographical 
areas. Other temporal variations in surface salinity exist at high lati- 
tudes where ice formation or melting alters the salinity profile. How- 
ever, the extent, doth vertically and horizontally, to which these 
varlations influence the open ocean are not known and can not be ade- 
quately modeled (Sverdrup et al., 1942). 

Allowing salinity to vary only latitudinally was dictated primarily 
by simplicity and econemy; however, the dependence in Y alone is justi- 
fiahle. The significant gradients in salinity are north-south and ap- 
preciable deviations from this pattern are found only at the boundaries 
where land runoff and boundary currents become important (Figure 1) 
(Riley and Skirrow, 1975). Since the model is primarily concerned with 
open ocean circulation, the assumption of longitudinal independence in 
Salinity seems reasonable. An anticipated result of the inclusion of 
Salinity into the model is the formation of an improved and more realis- 
fic subarctic gyre. This gyre could be created by the geostrophic flow 
induced by a north-south salinity gradient. 

The salinity cross section prescribed in the model was derived from 
actual measurements from research vessels. The series, Oceanic Observa- 


GiemseoLetne Pacitac: (Year) (Scripps Institution, 1960; 1963a and b; 
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1965a,b,c and d) compiled by Scripps Institution of Oceanography was 
the source of the measurements. Salinity profiles were selected to ap- 
proximate a cross section along 160°W longitude and at the model's grid 
point spacing of 2.03° latitude. Due to the sparseness of data in some 
areas, a common time of year could not be selected; however, 75% of the 
Measurements lie between late spring and early fall. Because the Pacific 
Ocean narrows and shoals above 60°N, data in this area are limited in 
both norizontal and vertical extent. To overcome the data insufficien- 
cles the profile at 60°N was duplicated for the grid points between 60°N 
and the model's boundary which is located at 65°N. 

Figure 2 shows the latitude-height cross section, taken from the 
above data sources, plotted and analyzed at the model's grid points. 
The salinity cross section (Figure 3) found in Intermediate Waters of the 
Pacific Ocean (Reid, 1965) was used to verify the representativeness of 
the compiled cross section. Major features such as the low salinity 
tongue penetrating equatorward from the high latitudes and the high 
Salinity central water mass are adequately represented. Areas of sharp 
Mom~ezontal Gradients at 35°N and 45°N are also reasonably depicted. 

Although the major features of the salinity cross section were well 
represented, a smoothing technique was employed to remove the small-scale 
irregularities and improve horizontal consistency. A repeated, three 


point weighted average given by 


Oo) 


t- 


was selected. Boundary values remained unchanged. This smoothing pro- 
cess was performed up to five times and the resulting cross sections 


analyzed. Four repetitions of this averaging scheme removed small-scale 


uD 





irregularities (less than two grid lengths) but retained the major fea- 
tures of the salinity structure (Figure 4). The central water mass re- 
Mains intact and the low salinity tongue is still identifiable. Hori- 
zontal gradients have been smoothed but maintain their location anda 
portion of their intensity. Since the cross section was composed of 
actual soundings, vertical consistency was assured and vertical smoothing 


was not considered necessary. 
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III. EQUATION OF STATE AND CONVECTIVE ADJUSTMENT 


The density of sea water as a function of temperature, salinity and 
pressure is usually computed using the Knudsen formula. Although the 
Knudsen formula is quite accurate, it is generally considered too inef- 
ficient for inclusion in numerical models. This inefficiency arises due 
to the wide range of temperature and salinity for which the Knudsen for- 
mula holds. Since the actual variation of temperature and salinity in 
the open ocean is known to be limited, several researcners have developed 
empirical formulas which are accurate and efficient but apply only to 
open ocean temperature and salinity ranges. Two such empirical formulas 
were considered for inclusion in the model. 

Bryan and Cox (1972) solved an overdetermined, simultaneous system 
of equations which related departures from the mean of temperature and 


salinity to the departure from the mean density at a given depth. That 
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The density is given by: 
x. OT + x. 6S + x (Sn) - + x (6s)° Poe eOhecu—mocr.c,Z )=Ox.)x 10° 
ie 2k ok Ak 5k aes Ok 


(3) 
The constant coefficients, x, were computed at each depth by Bryan and 
Cox uSing an iterative process. Acceptable accuracy was obtained with 
as few as three coefficients (n=3) and 50 temperature, salinity and den- 
sity combinations (m=50). A significant improvement in accuracy was ob- 
tained with nine coefficients (n=9, m=50). Figure 5 illustrates the 
accuracy of this method as compared to the Knudsen formula. The errors 
are always largest near the surface where the widest range of tempera- 
tures and salinities exist. 

Another empirical formula for the equation of state of sea water was 
developed by Friedrich and Levitus (1972). This formula is similar to 
the Bryan and Cox formula in that it fits a polynomial in temperature 
and salinity to computed values of density determined by the Knudsen for- 
mula; however, the coefficients are continuous, rather than discrete, 
Meet 1ons Of depth. 


A quadratic equation given by 


2 
= + + + 5 I ne Gaal 4 
Stes, 2.) aT cot oe Cy 2 S (4) 


is highly efficient but is not sufficiently accurate over the expected 
range of temperature and salinity. Further expansion of the polynomial 
co 

Z 


2 3 
= | Ee iy eC lL ST 5 
GS, las) c. + Cc, ~ c,s ~ c,T + C,ST C. _S (5) 
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ylelds acceptable accuracy while expansion to more terms does not im- 
prove accuracy appreciably. 

The coefficients are themselves quadratic polynomials in -Z (depth) 
which fit a smoothed curve to the pressure dependency of density. Due 
to the reduced range of temperature and salinity with depth, the coeffi- 
cients more accurately represent the oressure dependence at increased 
depth. Figure 6 shows the pressure dependence of the coefficients and 
the rms error of equation (5) with respect to the Knudsen formula. iIn- 
creased efficiency without a decrease in accuracy can be accomplished 
by reverting to equation (4) at depths of two kilometers and below. The 
values for the coefficients used in equations (4) and (5) are given in 
Tables I and II, respectively. Note that Z is in kilometers and defined 
positive downward. 

The Friedrich and Levitus equation of state was selected for inclu- 
sion in the model since computation of density at the model's irregular 
depths is facilitated by the continuous depth dependence. Another advan- 
tage over equation (3) of this formulation is that the equation of state 
would not have to be changed 1f the location of the levels in the model 
were ever changed. The explicit depth, temperature and salinity depen- 
dence of density as set forth by equation (5) aids in the remodeling of 
the convective adjustment process as discussed below. 

Without salinity the model had convectively adjusted whenever the 
profile was slightly positive, neutral or negative. The incorporation 
of salinity into the model enables a negative temperature profile to be 
convectively stable when coupled with a stabilizing salinity profile. To 
determine the relationship between depth, temperature, salinity and den- 


Sity, the total derivative was expanded into its components. 


JES; 





So Ee aS aS ae (6) 


soe do 
where Z increases upward. For stability, qz <9: the temperature pro- 


~ 


file must satisfy 


(i) Ceo et4). (5) 


The right hand side of equation (7) gives the temperature lapse rate 


which renders the density profile neutral (SS = QO). Terms (1) and (4) 
are known from the temperature and salinity profiles computed in the 
model. Terms (2), (3) and (5) are computed from the equation of state 
Memoeten (5) )2 For eiticiency, terms (2) and (3) can be simplified 
through scale analysis; however, all components of term (5) can be of 
Similar magnitude and must therefore be retained. The expansion of the 
terms in equation (7) and the scale analysis is detailed in Appendix A. 
iemede Inequality 1S v2olLated (as calculated using finite differencing 
between two levels), then the temperature profile is convectively ad- 
justed. The adjustment is accomplished by satisfying two constraints: 
1) the heat present in the two layers remains constant; and 
2) the new temperature profile equals the neutral temperature 
lapse rate, the right hand side of equation (7), plus a small 
additional stabilizing lapse rate. 


In this manner, cold water may overlie warm water provided the proper 


Seen Struceure Exists. 
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Vy hou GTS 


A "spin up" of 246 years of simulated time without salinity (Haney 
et al., 1978) was used as initial conditions. The salinity cross section 
and equation of state discussed above were inserted into the model and 
the integration continued for seven years. Streamlines and temperature 
contours at various depths for summer (mid-July) and winter (mid-January) 
of the model climatology, year 246 (without salinity) and the years 249 
through 253 (with salinity) were plotted, analyzed and compared. During 
the first seven year integration period, "neutral stability" was still 
defined as a e QO. That is, the stabilizing effect of salinity in the 
Geet ercrs2acion Of vertical convection (convective adjustment) was not 
included. This effect was added after year 253 (see below). It was 
handled in this manner in order to study the effects of one change at a 
Pemews lne model climatology, without salinity, Figures 7 through ll, 
depict the features previously discussed. In these and subsequent stream- 
line figures, the lengths of the small arrows are indicative of the rela- 
tive magnitude of the current at the grid point, while the long continuous 
lines are streamlines which help display the flow direction. In both the 
winter and summer the Equatorial Counter Current, North Equatorial Cur- 
rent, Kuroshio Current and North Pacific Current are well represented. 
At high latitudes a cyclonic turning is apparent but a closed Alaskan 
gyre is not formed. The temperature contours, Figures 12 through 16, 
illustrate the baroclinic zone formed by the North Pacific Current and 
the northward Gisplacement of this baroclinic zone by the Kuroshio Current. 


The contour interval is 3°C at all levels except at level 9 where the 
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Mieemval is 0.3°C. The North Pacific Central Water Mass is well defined 
iaenorezOntal and vertical extent by the 21°C isothern. 

After seven years of integration with salinity included, several 
changes in the streamline and temperature contours are visible. Compar- 
ing Figures 17 through 26 to their counterparts in the model climatology 
reveals a lessening of the cyclonic tendencies of the surface layers at 
high latitudes in the model with salinity included. Since this is 
directly opposed to the anticipated results, it merits further investi- 
gation. 

In the previous model with the equation of state represented by a 
linear function of temperature, the north-south gradients of temperature 
near 0°C, although small, generated significant geostrophic currents. 
Wind generated currents provided the cyclonic tendencies in the surface 
layers. The present equation of state is only weakly affected by the 
temperature fluctuations at low temperatures. The salinity gradients 
at high latitudes are very weak (recall, there is no horizontal salinity 
gradient from 60°N to 65°N) with the resulting current primarily induced 
By the wind fields (Figure 27). A strongly zonal forcing during winter 
generates meridional surface currents at high latitudes with the converse 
applying during summer. Below the surface, zonal currents resulting 
from the north-south temperature and salinity gradients are encountered 
but their magnitudes are small. 

Another feature of interest, which appears during winter (Figures 
17a through 21a) when salinity is included, is the southward directed 
current along the westem boundary near 45°N. This current 1s physically 
located in the region of the Oyashio Current. It appears only in winter, 


and since salinity is constant all year, atmospheric forcing is the 
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degical cause. The transport stream function (the vertically integrated 
mass transport is nondivergent and therefore represented by a stream 
function in this “rigid 1id" model) for winter as camnpared to summer 
(Figure 28) explains the seasonal dependence of this gyre. A strong cy- 
clonic stream function is located in this region in the winter while the 
value of the stream function in this area is small during summer. The 
salinity cross section (Figure 4) below 125 meters enhances the cyclonic 
shear due to a local salinity maximum coincident with the axis of the 
stream function minimum. The resulting streamlines depict a geostrophic 
gyre below 100 meters and a wind induced gyre above 100 meters. When the 
cyclonic atmospheric forcing diminishes in the summer, the upper levels 
revert to a zonal flow and the gyre disappears except at the lower levels. 
Penal ea change moted wich the inclusion of salinity is a rather large 
cyclonic vortex in the eastern portion of the domain at about 35°N. A 
time series of streamlines from July 248 through July 251 (Figures 29 
through 32) show that the vortex is oropagating westward at a rate of 
approximately 1.6 cm/sec. The vortex 1s present both winter and summer 
and penetrates from the surface to approximately 350 meters. The charac- 
teristics of this vortex point to an internal Rossby wave. The phase 


speed of an internal Rossby wave 1s given by 


ODS 


c= Ur (8) 


whe re U = speed of the mean current 


1, k and m are the wave numbers in tne x, y and z directions 
and 1 = 27/L 


£ = average value for the Coriolis parameter 
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Lee) 
il 


change in f with latitude 


oT, 
N= (ag aa the Brunt-Vaisala frequency 


This equation is based on a constant mean flow, U, and a Brunt-Vaisala 
frequency, N, which is qualitatively applicable to the region in question. 
Values for the above equation were determined from the streamline plots, 


temperature profiles and velocity fields. The following values were 


selected: 
; ~5 -1 
f mee, Sin = 20) sin 3555 = 6.4 x 10 sec 
a A = -13 a ee 
= = 36 f 35 io x eo sec cm 
ieloe on 
Ll = 27/1000 km 
Kee=a2n73sc00 km 
m= 27/800 m 
2 4 4 -l 


No = 2x10. x 980 x 8/2250 cm = 7x10 sec 

Ue= 1.5 cm/sec 
Inserting these values into equation (8) yields a value for the phase 
speed of -1.5 cm/sec which compares favorably with the observed speed. 
At this speed the wave would propagate across the domain of the model in 
12 years. 

The effects of the Rossby wave are well illustrated by the changes 
in the temperature contours (Figures 33 through 36). With the approach 
of the wave from the east, the baroclinic zone is displaced southward. 
fame scnie Daroclinic zone nearly returns to Lts original location after 
the wave passes, the temperature structure of the high latitudes is al- 
tered for a more extensive period. The high latitude lsotherms are dis- 
placed southward and eastward; however the streamlines indicate flow from 
the south which should aid in returning the isotherms to their previous 


Genel guration. 
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The Rossby wave is therefore considered to be a transient feature 
due to the sudden insertion of salinity into the model and the salinity's 
dynamical effects on the currents. As the wave slowly propagates across 
the basin, its amplitude decreases due to the lateral eddy diffusion in 
the model. The e-folding time scale for damping is of the order of 

2 2 -l : 
((k° + 1°) AJ , where ae 1s the lateral eddy diffusion for heat. With 

7. y) . 

A fre emi /sec, the time scale for damping is of the order of 5-10 years. 
Thus it appears that another 5-10 years of model integration will perhaps 
be necessary in order to allow the large-scale GEOSe poe eucueren ts -cO 
completely adjust to the new dynamical constraints imposed by the salinity 
field. This is not considered a serious problem. 

With the dynamical effects of the salinity reasonably well defined, 
the convective adjustment parameterization discussed in the prior section 
was inserted into the model. Within 40 days of simulated time tne tempera- 
tures exceeded the Fr6.2 format. The individual terms of equation (7) were 
rechecked to insure their magnitudes and signs were correct. All terms 
produced appropriate values for realistic salinities, depths and tempera- 
tures; however, in the presence of strong surface cooling at high latitudes 
during the winter of year 253-254, unrealistic sea temperatures were 

| J0 : or ; 
generated. Temperatures below -3°C caused 3p tO become positive which 
oT ae 
forced the neutral value of az to become large and positive. The pcara- 
. 

meterization of convection permitted temperatures below -3°C; however, 
the equation of state does not hold for temperatures below -3°C and dces 
not model the change of phase from water to ice which occurs at subzero 
temperatures. 


While distilled water reaches maximum density before freezing, sea 


water with salinities greater than 25% exhibit increasing density with 
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decreasing temperature throughout the cooling process; therefore, super- 
cooled water continues to be convectively mixed (Sverdrup et al., 1942). 
Should the water column become supercooled through its entire depth, an 
ice cover would form insulating the column against further heat loss. 

The convective adjustment parameterization was altered to prevent 
the instabilities previously encountered. Water cooled below -2°C was 
convectively mixed so as to conserve heat but maintain minimum tempera- 
Mmiees ae © above -2°C. Should the entire column cool to -2°C no further 
cooling would occur. With these changes inserted, the model was run 
through a simulated winter. Surface temperatures rapidly dropped to 
-2°C at the higher latitudes and were consistently cooler than the under- 
lying layers, the stratification being provided by the imposed salinity. 
Because low temperatures nave little effect on density, changes in the 
Circulation voatterns were not noted. The inclusion of a convective ad- 
justment orocess dependent on depth and salinity as well as temperature 
allows cold anomalies to exist and more accurately models the physical 


(convective) processes observed in the high latitude ocean during winter. 





Vie CONCLUSIONS 


The inclusion of a prescribed salinity cross section into the model 
produced limited changes in the large-scale ocean circulation patterns. 
The most significant change was the appearance of a cyclonic gyre in the 
northwest portion of the domain. While its seasonal fluctuation at the 
surface is due to atmospheric forcing, the local salinity maximum below 
G@Z5 meters induced cyclonic shear at this latitude. No cyclonic gyre 
was formed in the region of the Gulf of Alaska presumably due to the 
reduced magnitude of the cyclonic atmospheric forcing in the northeast 
portion of the domain. The lack of bottom topography and of irregular 
boundaries may also impede the formation of an Alaskan gyre since all cy- 
mone VOrcicity, which 1s put into that region by the wind, is free to 
propagate westward all the way to the western boundary. 

The propagation of an internal Rossby wave across the mid-latitudes 
altered circulation and isotherm patterns; however, these changes were 
Of a nonpermanent nature. The effect of dissipative forces were apparent 
after six years. Based cn the orescribed value of the lateral eddy dif- 
#usion coefficient of Seo. ee eee the Rossby wave will be dissipated 
in approximately 10 years. 

The convective adjustment parameterization allows negative tempera- 
Pimewproriles to exist when balanced by a stabilizing salinity profile. 
To prevent supercooling of the surface layers, a simplistic mixing and 
ice formation parameterization was required. Although reasonable results 
were obtained, a more thorough investigation of the process of vertical 


convective mixing at freezing temperatures is recommended. 
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Several problems were encountered because the model extends to 65°N. 
Geographically, the Pacific Ocean above 60°N can not be realistically 
represented by a basin 4 kilometers deep with a width of 90° longitude. 
Salinity data above 60°N are sparse and highly variable since land runoff 
and ice formation or melting significantly alters the salinity profile 
there. For these reasons, 60°N would be a more appropriate northern 
boundary for the model. 

While the inclusion of salinity in the model did not generate all 
the desired results, it has improved the model in specific areas. The 
inclusion of salinity and the improved convective adjustment process 
are logical progressions in the formulation of realistic ocean models 


which can be effectively used in the study of large-scale ocean dynamics. 
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APPENDIX A 
DETAILS OF THE CONVECLDIVE 
ADJUSTMENT PARAMETERIZATION 
The empirical equation of state formulated by Friedrich and Levitus 


(1972) is given by: 


2 3 2 
Os,t,2) = C. + Cc + + ~ ST + Me ee oC 
( ) 1 ot c,s c,T C. C. 4ST 


The values for the coefficients, which are polynomials in -Z (depth), 
are given in Table I. Note that Z is in kilometers and defined as posi- 
tive downward. Changing Z tc positive upward and taking the partial 


derivative of sigma with respect to Z yields: 


dC dC gC dC dC dC 2 8 oC_ 5 


5 
= aan — —- —— —- ——a_—— — li <—e ; Uh ar as am 5 ST 
32Z ja) Se Se eee oy > ee i 
Inserting the values for the coefficients and performing a scale analysis 


on the equation yields: 


dO 
pee gels  -1Z + 203639 — -O0lezT + .09865 - .0002SZ 
~.0006T* - .0002ST + 
Bs = cao 
See 215) (10 gm/cm km) 


=, 


Bey ieis x 10° (0. gm/em> cm). 


Likewise, the differentiation and scale analysis for sigma with respect 


Goms 1S Given by: 


Z 
aed 
c, +n Cot Cot 
2 
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8056 (10° gm/em”s) ; 
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however, the expression for the change in sigma with respect to T can 


not be simplified to a single value through scale analysis. Retaining 


all terms which may be significant over the expected range of tempera- 


ture, Salinity and depth, the effect of temperature on density is given 


iy: 


OT 


or 


where 


J v 


2 


2 
+ 2c,T oF C.S ao 3C,7 + 2ClsT 


2 


Mes + .Oe6a4+ Ol0eZ = .01527T — .0013Z2T - .003S = -8002Zs 


+ 00018" Tt OOO LST. 


stable lapse rate is given by: 


- apes Saree : ae. 
WSO -8056 =F) /(K, + KT + Ks - .QOOO1LT (T+S) ) 


2 
GS = .496)+e 02052 + .0008Z 
Ry = - 70152). 6.001G4 
RK, = =-,.003 — .00022 
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TABES 1. 
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Coefficients for the Friedrich and Levitus equation 

Ci stace meguactonm(s)e) GC. = 4,0 4 5.2 + c.Z* where 
are een: a 

Zis in kilometers and positive downward. 


a Igy eS 
a a Jib 
ay -0 =2 
-7.2169x10 +5.1215xL0 -5 .012x10 
+4.9762x10_ -3.6349x10 ~~ +7.853x10_, 
+8.0560x10_ -8.5540x10_, +1.070x10 
=7. 5911x102 +6.4295x10_/ -1.397x10_> 
-3.0063x10_ +1.9365x10 -3.899x10_¢ 
+3.5187x10_, ~3.9740x10_ -5 .695x10_ 
+3.7297x10 -2 .8108xL0 +1.147x10 


Coefficients for the Friedrich and Levitus equation 


of state, equation (4). C. =a, +b.Z + c,Z° where 
25S (Aneel omerers and positive downward. 
an bg or 
a af af 
=> -0 =2 
-9.2163x10_, +5.1140x10_ -4.692x10_7 
+4.3314x10_) =o} 2lah=jop.dliels +7.689x10_ 
+8.0640x10_, ~8.S826x10_), +1.433x10_- 
=62 2725-0) +5.1351x10_, -1.246x10_ 
=2-/762x00 eo 2 0 =3.985x10 


Si: 





“(tP6t '°Te 22 dnzpz9as) ues0Q OTF TOeRd By JO SsANO\UOD 


ie 
a 
Ad / ee 


ve 


Ajtutypes soezans 





yp Shae cay 


yan rafaid oaso yonby 
WNSOPOIUOYT FAPODD 





SZ 











tes 
-—: ca 
ee ee yaa 
Si 
\ 


Be 


- a 


Depth Gn) 


. Depth (ec 


‘ 
el ee eee ey, 


ou = 





er 
—— eee SE ee A ES TS 


Figure 2. Salinity cross section along 160°W with no horizontal 
smoothing. 
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@igure 4. Salinity cross section along 160°W smoothed four times. 
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Figure 5. Average percentage error in sigma for the Bryan and Cox 
equation of state where the Knudsen formula is taken as 
standard (Bryan and Cox, 1972). 
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Figure 6. Coefficients of the Friedrich and Levitus equation of state 
as a function of depth and rms error of this equation with 
respect to the Knudsen formula (Friedrich and Levitus, 1972). 
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Figure 7. Streamlines and velocity vectors for surface layer, year 246 
(without salinity) winter (top) and summer (bottom). 
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Figure 8. Level 3 streamlines and velocity vectors 
year 246 (without salinity) winter Crop) sand 
summer (bottom) 
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Figure 9 Level 5 streamlines and velocity vectors 
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year 246 (without salinity) winter Ctop) and 
summer (bottom) 
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Figure 10. Level 7 streamlines and velocity vectors 
year 246 (without salinity) winter (top). sand 
summer (bottom) 
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Figure 11. Level 9 streamlines and velocity vectors 
year 246 (without salinity) winter (top) and 
summer (bottom) 
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Figure 12. Temperature contours for surface layer, year 246 (without 
salinity) winter (top) and summer (bottam). 
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EPaspresia. Level 3 temperature contours, year 246 
(without salinity) winter Ctop) and summer (bottom) 
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Figure 14. Level 5 temperature contours, year 246 
(without salinity) winter (top) and summer (bottom) 
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Figure 15. Level 7 temperature contours, year 246 
(without salinity) winter (top) and summer (bottom) 
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Figure 16. Level 2 temperature contours, year 246 
(without salinity) winter (top) and summer (bottom) 
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Figure 17. Surface streamlines and velocity vectors 
year 253 (with salinity) winter (top) and 
summer Cbottom) 
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Figure 18. Level 3 streamlines and velocity vectors 
year 253 (with salinity) winter (top) and 
summer (bottom) 
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Figure 19. Level 5 streamlines and velocity vectors 
year 253 (with salinity) winter (top) and 
summer (bottom) 
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Figure 20. Level 7 streamlines and velocity vectors 
year 253 Cwith salinity) winter (top) and 
summer (Cbottom) 
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Figure 22. Surface temperature contours, year 253 
Groene salinity) winter C€too) and summer (bottom) 


ay 





3 0 
145E 155 165 175 180 175 165 155 145 135 125W 











2 ee 


3 | ae O 
145E 155 165 <175180175 165 155 1453s 135 125W.| 





a A I I Ry re UP er 


igure 23. Level 3 temperature contours, year 253 
(with salinity) winter (top) and summer (bottom) 
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Figure 24. Level 5 temperature contours, year 253 
(with salinity) winter (top) and summer (bottom) 
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Pee. bevel) Stemperature comeours, year 253 
(with salinity) winter (top) and summer (bottom) 
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Figure 26. Level 9 temperature contours, year 253 
(with salinity) winter (top) and summer (bottom) 
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Figure 27. Climatological surface wind stress 
patterns winter (top) and summer (bottom) 
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Pigure 28 Transport st I l 
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Figure 37. Climatological atmospheric quantities used in the 


model. (a) Solar radiation at the top of the 
(Wm ~); (b) fractional cloud cover (tenths); 
temperature (°C); (d) surface vapor pressure 
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